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We present a new dynamical approach for measuring the temperature of a Hamiltonian dynamical system in the micro canonical 
ensemble of thermodynamics. We show that under the hypothesis of ergodicity the temperature can be computed as a time- 
average of the functional, V • (Vi//1| VHj| 2 ), on the energy-surface. Our method not only yields an efficient computational 
approach for determining the temperature it also provides an intrinsic link between dynamical systems theory and the statistical 
mechanics of Hamiltonian systems. 
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Two crucial assumptions allow one to pass from the 
dynamical to the thermodynamical description of a clas¬ 
sical Hamiltonian system. First, the system should ex¬ 
hibit some sort of ergodicity so that time-averages can be 
replaced with equivalent space-averages over the micro- 
canonical ensemble, i.e. over the corresponding energy 
surface (possibly with other constraints) and with respect 
to the invariant Liouville measure. Second, it should 
make sense to consider a large number of (weakly) cou¬ 
pled sub-systems and look at the statistical behavior of 
an individual fixed sub-system. As the size of the ambient 
system tends to infinity one obtains the thermodynamic 
limit in which the micro-canonical ensemble involving the 
sub-system is assumed to behave like the canonical or the 
Gibbs ensemble where quantities like temperature and 
pressure become operationally defined. In a few special 
cases it has been possible to prove the equivalence of the 
two ensembles but it is believed that under quite general 
assumptions the two ensembles indeed give equivalent de¬ 
scriptions when systems are large. We refer to Ruelle Q 
for a more complete treatment on the equivalence of en¬ 
sembles, to Lebowitz et al. || for an illustrative example 
of some differences between the ensembles and to Landau 
and Lifshitz j2| for a general introduction to the subject. 
For practical calculations using the micro canonical en¬ 
semble see e.g. Evans and Morriss |jj. 

In the thermodynamic description the temperature 
plays a prominent role and the problem of measuring it 
using the dynamical evolution of the system dates back to 
the early days of thermodynamics. Already Maxwell re¬ 
alized that when the Hamiltonian has the special form 
H(p,q,...) = p 2 /2 + K(q, ■ ■ •), thus singling out the 
square of a momenta, the canonical ensemble average of 
p 2 happens to be the temperature. Thus assuming equiv¬ 
alence of ensembles and ergodicity it suffices to measure 
the time-average of p 2 during the dynamical evolution in 
order to get the temperature. 

Here we shall pursue a different approach in which we 
focus on the global geometric structure of the energy sur¬ 
face We shall assume ergodicity only and show that the 
temperature as defined in the micro-canonical ensemble 


can be obtained as a purely dynamical time-average of a 
function which is related to the curvatures of the energy 
surface. 

In this presentation we shall for clarity and without 
loss of generality restrict ourselves to a Euclidean phase 
spaced, Q = R 2d (d > 1), and a Hamiltonian function, 
H : fl —*■ R, which preserves the Liouville measure (here 
Lebesgue measure in) and such that for all energies con¬ 
sidered the set H(£) < E is bounded and has a smooth 
and connected surface, A(E) = {H(£) = E}. Further¬ 
more, we shall consider the case when the energy itself, 
H = E, is the only preserved first integral. We refer to 
Q for a generalization of this approach. 

For the micro-canonical ensemble of classical thermo¬ 
dynamics one defines the exponential of the entropy to 
be the (canonically invariant) weighted area of an energy 
surface : 

eStm s /’»«> nmo-m = <u 

and the temperature is obtained through differentiation 


T{E) dE' w 

In practise a calculation of the entropy function is im¬ 
possible when the number of degrees of freedom in the 
system is large. However, we have the following : 

Theorem : Suppose that (fi, H) is ergodic with respect to 
the Liouville measure m\A(E), restricted to a non-singular 
energy surface, A(E). Then for m\A(E)-almost every 
a;(0) £ A(E), one has 

—J— = lim j f dr (3) 

T(E) t—°o t J o 


*The Euclidean metric is not necessary but provides a con¬ 
venient way of expressing the dynamical average in terms of 
gradients of functions and divergences of vector fields. 
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where x(r) is the trajectory of ir(0) under the Hamilto¬ 
nian flow and the observable $ equals V • (ViL/1| VHj| 2 ), 
which can also be written as : 


is the Hamiltonian flow and g j) is the normalized energy 
gradient flow. Since = (VH(x),r](x)) = 1 

we have 


AH n D 2 H{\7H,VH) 

ww ww 


Remarks : 

1. The measure and the energy are canonical invari¬ 
ants. Hence, from equations 0 and ||), so are the 
temperature and the time-average, equation !), of 
$ but not $ itself. The theorem shows also that the 
temperature unlike the entropy is an intrinsic dy¬ 
namical feature of the system (it can be computed 
as a dynamical average). 

2. The theorem can quite easily, at the cost of making 
the notation more heavy, be generalized to arbi¬ 
trary symplectic manifolds without a Riemannian 
structure (cf. ||). 

3. It is of interest to see how the conjugated thermo¬ 
dynamic variables of other preserved first-integrals 
(e.g. total angular momentum) can be incorpo¬ 
rated in the above dynamical formalism. It is for 
example well-known that the pressure can be ob¬ 
tained as a time-average of the momentum trans¬ 
fer per surface area element in the configuration 
space. It is also of interest to derive general 
(thermo-)dynamical relations for such quantities in 
the micro canonical ensemble (cf. Q). 

4. There is a formal resemblance of <f> to the Gauss 
curvature 1/R = V • (S7H/\\ VHj|), where for the 
latter quantity the norm of the gradient is not 
squared. In particular, when the number of de¬ 
grees of freedom gets large (and the last term in 
the formula ([|) becomes relatively small) the rela¬ 
tionship becomes closer. We do not know of any 
deeper reasons for this connection. 

The idea behind the above Theorem is to cast the en¬ 
ergy dependence of the geometry, i.e. the dependence on 
A(E ) in equations 0 and (||), into an algebraic depen¬ 
dence for a fixed geometry. This is achieved through the 
introduction of an auxiliary vector field (to a large extend 
arbitrary) and it will enable us to take the derivative with 
respect to the energy inside the integral. 

We consider in the following two vector fields on H : 

„ = JVff and (5) 

Here J is a symplectic matrix, e.g. ^ i o ) ’ an< ^ ^ 

is the usual Euclidean product in fl. For a vector field 
I on SI we let g x denote the associated flow. Thus gl 


h (9^(x)) = s + H{x). ( 6 ) 

Assuming that no singular points are encountered along 
the gradient flow, i.e. that S7H 0 everywhere on the 
energy surface (which is a generic condition), the map 


9q : A{E) —> A{E + s) (7) 

is a diffeomorphism between energy surfaces. 

Both the Liouville measure rn and the energy are pre¬ 
served under the Hamiltonian flow, i.e. m(g t v dx) = 
m{dx) and H(gl(x)) = H(x). It follows that by con¬ 
sidering an infinitesimal flow-box (under the energy gra¬ 
dient flow) g^'^dA and by measuring the volume of such 
a flow-box we may define a measure g on the energy sur¬ 
face A{E ) which is also invariant under g f v , 

g(dA) = lim -m{q^^dA\. (8) 

e—> 0 + e ' 


For this measure, we have = J A ^ E ) g(dA). The flow 
g^ commutes with itself and hence it maps the flow-box 

g\®'^dA(E) into the flow-box g^'^dA(E+s) (by equation 
[?]). But then the ratios between the weighted areas of in¬ 
finitesimal surface elements just becomes a ratio between 
volume elements which is nothing but the Jacobian of the 
mapping : 


g{dA{E + s)) 
KdA(E)) 


do s 

= det(—--) = l + sV'ij + 0(s 2 ) 
ox 


(to order s). Defining $ = V ■ g (i.e. our equation Q) we 
get the relation : 


[ g(dA) = f g(dA)(l + + ...), 

JA(E+s) J A(E) 

and finally, using the definition of temperature 
obtain : 

1 Sa(e) 9 {dA)<& 

T{E) = f A(E)f i(dA) ■ 

The interesting observation is now that if the Hamilto¬ 
nian flow is ergodic with respect to the Liouville measure 
restricted to the energy surface, g, then by Birkhoff’s 
theorem the right hand side equals the time-average (and 
thus equation [|) of $ for /r-almost every (equivalent to 
m\A(E)- a - e -) trajectory under the Hamiltonian flow, thus 
proving our theorem. 


(9) 

d) we 
(10) 
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An illustrative, though trivial!!], example is given if 
we take N harmonic oscillators with the Hamiltonian 
H = YjiLi + Pi)- I ts gradient \7H = (q,p) is non¬ 
vanishing at all non-zero energies and a straight forward 
calculation yields <& = (iV — 1 )/E = 1/T(E). In this par¬ 
ticular case we can directly read off the micro-canonical 
result, E = (N — l)T. The fixation of the energy ‘freezes’ 
one degree of freedom as is seen by comparing with the 
canonical result, (E) = NT. 

In conclusion we have described a dynamical way of 
measuring the temperature in the micro canonical en¬ 
semble of thermodynamics. It is rigorously shown to 
work for a Hamiltonian system at energies where the en¬ 
ergy surface is regular and the flow is ergodic. The ap¬ 
proach is purely dynamical and involves calculating the 
time-average of the functional, V • (VH/\\S7H || 2 ), on the 
energy surface. As the derivation shows the dynamical 
temperature is really intrinsic to the dynamical system in 
contrast to the entropy itself. We note that in a classical 
theory the entropy is in any case only defined up to a 
constant. We believe that the results and the underlying 
ideas presented here may be of practical importance in 
numerical simulations and helpful to physicists who wish 
to understand the foundations of classical thermodynam¬ 
ics. 

I am grateful to Predrag Cvitanovic, Robert 
MacKay, David Ruelle as well as to many referees for 
(thermo-)dynamic discussions. 


[1] D.E. Evans and G.P. Morriss, Statistical Mechanics of non- 
Equilibrium Liquids, Acad. Press (1990). 

[2] L.D. Landau and E.M. Lifshitz, Statistical Physics, Part 1, 
Oxford, Pergamon (1980). 

[3] J.L. Lebowitz, J.K Percus and L. Verlet, Ensemble Depen¬ 
dence of Fluctuations with Application to Machine Compu¬ 
tations, Phys. Rev. 153 , 250 (1967). 

[4] D. Ruelle, Statistical mechanics: Rigorous Results, 
Addison-Wesley (1989). 

[5] H.H. Rugh, A Dynamical Approach to Thermodynamics in 
the p-Canonical Ensemble (in preparation). 


t Strictly speaking this example does not meet the conditions 
of our Theorem as it is not ergodic. 
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